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Introduction 


The  following  discussion  repeats  that  reported  in  Ref.  2. 

Tne  wake  of  a  helicopter  rotor  consists  of  &  relatively  concentrated 
tip  vortex  generated  by  the  rapid  decrease  in  bound  circulation  at 
tne  tip  5f  the  rotor  blade  and  distributed  vortex  sheets  generated 
by  spanwlse  variation  of  the  bound  circulation  over  the  inboard  per— 
tion  of  the  blade  (the  inboard  trailing  wake)  and  by  azimuthal  vari¬ 
ations  of  the  bound  vorticity  (the  shed  wake).  The  tip  vortex,  be¬ 
cause  of  its  concentration,  is  best  represented  by  a  vortex  line, 
with  a  finite  vortex  core  where  viscous  effects  dominate,  and  is 
responsible  for  the  sharp  peaks  in  the  rotor  airloads  distribution. 
Tnus  accurate  geometry  is  much  more  important  for  the  tip  vortex 
than  for  the  rest  of  the  wake  and  effort  has  been  concentrated  on 
the  accurate  and  efficient  computation  of  the  tip  vortex  geometry 

The  computation  of  the  tip  vortex  geometry  involves  the  in¬ 
tegration  over  time  of  the  induced  velocity  at  various  points  on 
the  tip  vortex.  Since  the  computed  tip  vortex  geometry  changes  as 
the  computation  proceeds  (converging  on  the  actual  geometry),  the 
induced  velocities  must  be  recomputed  (updated)  periodically  to 
account  for  these  changes.  Since  moot  of  the  computational  effort 
goes  into  computing  induced  velocities,  finding  techniques  to 
minimize  the  amount  of  updating  required  i3  very  important. 

Any  point  P  on  the  tip  vortex  can  be  Identified  by  its 
azimuth  and  its  age  (time  elapsed  since  point  P  was  trailec.  from 
the  rotor)  at  any  instant.  Tip  vortex  geometry  computations  gen¬ 
erally  proceed  by  computing  the  geometry  of  every  azimuth  at  age 
AJ>  (typically  n  15°),  then  at  2A<j>,  3Ap,  and  so  on  up  to  Jjtt  or 
.  This  computation  is  repeated  iteratively  until  satisfactory 
convergence  is  achieved.  In  Reference  1.,  the  induced  velocities 
used  tc  compute  the  geometry  were  updated  every  nA<J>  in  ago 
(typically  n  *  6)  at  every  azimuth.  At  every  age  at  which  In¬ 
duced  velocities  were  updated  the  geometry  for  ages  from  A<i-  to  & 
and  at  all  azimuths  was  recomputed  to  reflect  the  updated  Induced 
velocities.  A  technique  for  updating  the  induced  velocity  ana 


1 


geometry  at  each  azimuth  independent  of  the  other  azimuths  based 
cn  the  amount  of  distortion  accumulated  was  developed  under  con¬ 
tract  N00019--67-C-0236  and  will  be  described  in  Reference  2.  This 
net hod  of  independent  updating  did  not  result  in  the  desired  re¬ 
duction  in  computational  effort  and  was  abandoned  for  work  on  more 
promising  techniques  under  contract  N00019-68-C-0150,  as  described 
below . 

Now  Work 

The  following  discusses  work  beyond  that  reported  in  Reference 
The  basic  idea  was  to  use  an  updating  technique  similar  to  that  of 
Reference  1  where  the  induced  velocities  and  geometry  at  all  azimuths 
were  updated  every  nA$  in  age,  but  to  update  only  the  induced  veloci¬ 
ties  generated  by  part  of  the  wake  (the  near  wake)  and  leave  the  in¬ 
duced  velocities  generated  by  the  rest  of  the  wake  (the  far  wake) 
fixed  during  each  iteration.  The  far  wake  must  be  chosen  to  be  so 
far  away  from  the  point  P  at  which  induced  velocities  are  computed 
that  normal  changes  in  the  computed  wake  geometry  are  negligible  by 
comparison.  On  the  other  hand,  the  near  wake  should  be  minimized  to 
reduce  trie  computational  effort  required. 

The  induced  velocity  generated  by  the  entire  wake,  both  near 
and  far,  are  computed  once  at  the  beginning  of  each  iteration  at 
every  point  P  on  the  tip  vortex  at  which  the  wake  geometry  is  to  be 
computed.  Thus  the  ideal  way  to  divide  the  wake  is  to  define  the 
near  wake  relative  to  any  point  P  to  include  any  element  of  the  wake 
model  which  generates  an  induced  velocity  at  point  P  greater  than  V  , 
an  input  parameter  to  be  determined  by  experiment.  Since  there  can 
La  up  to  182?  points  P  and  up  to  1750  elements  of  the  wake  model 
which  generate  Induced  velocity  at  any  pcint  P  about  3  x  10^  storage 
locations  would  be  needed  to  store  such  a  map  of  the  near  and  far 
wakes.  Since  tne  available  computer  facility  has  about  4  x  10^  words 
of  core  storage,  this  method  of  defining  the  near  wake  was  temporarily 
abandoned  and  various  attempts  were  made  to  develop  a  simple  a  priori 
definition  cf  the  near  wake.  These  included  defining  all  of  the  tip 


vortex  to  be  the  near  wake  or  defining  all  of  the  wake  within 
(typically  n/2  or  tt)  in  azimuth  of  any  point  P  and  on  the  3ame 
tip  vortex  as  that  point  P  to  be  the  near  wake.  None  of  these 
simple  definitions  of  the  near  wake  gave  acceptable  results  compared 
to  the  normal  computation  without  near  and  far  wakes. 

At  this  point  in  time,  a  review  copy  of  the  paper  which 
eventually  became  Reference  3  was  received  and  it  reported  the 
very  successful  application  of  the  ideal  definition  of  the  near  wake 
described  above.  Since  the  core  storage  needed  to  use  this  defini¬ 
tion  of  the  near  wake  was  not  available  at  M.I.T.,  it  was  decided 
to  attempt  to  develop  an  approximation  to  the  ideal  definition  of 
the  near  v.'ake  which  would  minimize  storage  requirements.  The  first 
approximation  involved  gathering  all  of  the  wake  elements  at  each 
azimuth  and  age  combination  on  the  wake  of  each  blade  into  a  unit. 

Each  unit  included  an  element  of  the  tip  vortex,  an  element  of  the 
inboard  trailing  vortex  sheet,  and  an  element  or  elements  of  the 
shed  wake.  The  near  wake  relative  to  each  point  P  on  the  tip  vortex 
was  redefined  to  include  every  unit  which  generated  a  total  induced 
velocity  at  point  P  greater  than  V^. 

If  unlimited  core  storage  were  available,  the  easiest  thing 
to  do  would  be  to  store  the  induced  velocity  generated  by  each 
unit  at  each  point  P  on  the  tip  vortex  in  a  separate  storage  loca¬ 
tion.  The  computation  of  the  updated  induced  velocity  at  each  point  P 
on  the  tip  vortex  would  then  have  proceeded  by  examining  induced  ve¬ 
locity  at  point  P  generated  by  each  unit.  If  the  induced  velocity 
was  greater  than  Vm  it  would  be-  recomputed,  otherwise  the  old  induced 
velocity  would  simply  be  added  on  to  the  total.  Since  the  core 
storage  required  for  this  was  not  available,  the  total  induced  velocity 
at  each  point  P  (each  different  azimuth  and  age)  on  the  tip  vortex 
generated  by  the  far  wake  (consisting  of  all  units  generating  an  in¬ 
duced  velocity  less  than  at  point  P)  was  computed  once  at  the  be¬ 
ginning  of  each  iteration  and  stored. 

A  map  of  the  near  and  far  wakes  was  also  needed  since  this 
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information  was  no  longer  contained  in  the  induced  velocity  in¬ 
formation.  Instead  of  storing  a  different  map  of  the  near  and 
far  wakes  for  each  different  point  P  (each  different  azimuth  and 
age  combination)  on  the  tip  vortex,  a  different  map  was  stored 
for  each  different  azimuth  on  the  wake  of  each  blade.  This  saved 
storage  since  there  are  fewer  blades  than  different  ages.  The 
maps  were  stored  in  terms  of  the  time  elapsed  between  *he  instant 
when  the  element  of  the  tip  vortex  at  point  P  was  trailed  by  the 
rotor  and  the  Instant  when  the  unit  of  wake  which  is  to  be  classi¬ 
fied  as  near  or  far  wake  relative  to  point  P  was  shed  by  the  rotor. 

To  save  storage  only  the  elapsed  times  for  the  transitions  between 
near  and  far  wake  were  stored  Instead  of  storing  the  elapsed  times 
for  all  of  the  units  in  the  near  wake. 

The  net  result  of  the  approximations  described  above  was  to 

reduce  the  core  storage  required  to  use  the  near  find  far  wake  model 
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from  about  3  x  10  words  to  about  1  x  10  words.  This  made  it 
practical  to  develop  a  new  computer  program  using  the  near  and  far 
wake  approximation.  The  large  reduction  in  the  amount  of  computation 
involved  in  updating,  due  to  the  use  of  the  near  and  far  wake  approxi¬ 
mation,  allowed  the  basic  updating  scheme  of  Reference  1  to  be  modi¬ 
fied  for  Greater  accuracy  without  resulting  in  prohibitive  amounts  of 
computation.  In  Reference  1,  the  updating  of  the  induced  velocity 
and  the  wake  geometry  at  any  age  <5  and  azimuth  ^  was  done  by  re¬ 
computing  the  induced  velocity  generated  by  the  entire  wake  at  age 
zero  and  azimuth  ^  and  then  adding  on  the  induced  velocity  generated 
at  azimuth  4/  by  each  additional  unit  of  wake  that  was  created  as  the 
age  was  increased  in  steps  of  up  to  6 .  This  gave  a  tine  history 
of  the  induced  velocity  at  azimuth  ^  for  all  ages  between  zero  and 
5  in  A<J>  steps  which  was  integrated  to  get  the  corresponding  time 
history  of  wake  geometry.  To  improve  the  accuracy  of  the  computation, 
this  updating  technique  was  modified  by  recomputing  the  induced  ve¬ 
locity  generated  by  the  entire  wake  at  azimuth  and  age  0,  then 
adding  on  the  induced  velocity  generated  at  azimuth  \p  by  each 
additional  unit  of  wake  created  as  the  age  was  increased  in  steps  of 
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A<j>  up  to  (n-l)A$,  then  recomputing  the  induced  velocity  generated 
by  the  entire  wake  at  azimuth  ^  and  age  n&$,  and  so  on  up  to  age  6. 
Despite  the  large  amount  of  additional  recomputation  of  the  Induced 
velocities  Involved  in  this  improved  updating  technique,  the  new 
wake  geometry  computer  program  developed  using  it  plus  the  near  and 
far  wake  approximation  proved  to  be  substantially  faster  than  the 
program  described  in  Reference  1. 

The  new  wake  geometry  program  was  tested  using  various  values 
V_  and  n.  Preliminary  results  indicate  V  ■  0.03  and  n  *  8  give 
good  results  but  more  testing  is  necessary.  When  compared  with  the 
wake  geometry  program  described  in  Reference  1,  the  new  wake  geo¬ 
metry  program  was  not  only  faster  but  gave  somewhat  different  re¬ 
sults.  In  theory,  the  results  from  the  new  program  should  be  more 
accurate  than  those  obtained  in  Reference  1,  but  uncertainties  in 
the  experimental  wake  geometry  data  from  smoke  tests  prevented  con¬ 
firmation  of  this.  Work  is  continuing  under  contract  N00019-69-C- 
0219  to  refine  the  experimental  data  to  provide  a  reliable  check  on 
the  theory  and  to  refine  the  nev  wake  geometry  program.  This  work, 
including  a  more  detailed  description  of  the  new  wake  geometry  pro¬ 
gram,  will  be  reported  in  Reference  2. 

In  parallel  with  the  development  of  the  new  wake  geometry 
program,  work  continued  on  the  U3e  of  vortex  sheets  to  represent 
the  Inboard  trailing  wake  and  the  shed  wake  in  the  wake  model.  This 
is  desirable  because,  as  discussed  in  the  Introduction,  the  inboard 
trailing  wake  and  the  shed  wake  are  distributed  vortex  sheets  in 
the  real  world.  When  vortex  lines  are  used  to  represent  these 
distributed  vortex  sheets,  artificial  peaks  are  introduced  in  the 
induced  velocity  distribution  and  it  is  a  matter  of  chance  whether 
or  not  a  point  P  on  the  tip  vortex  at  which  induced  velocities  rre 
computed  happens  to  lie  on  the  positive  or  on  the  negative  induced 
velocity  side  of  such  a  peak.  The  use  of  vortex  sheets  eliminates 
this  problem. 

A  closed  form  integration  of  the  Biot  Savart  relation  for 
induced  velocities  over  a  rectangular,  planar  vortex  sheet  has  been 
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developed  similar  to  that  developed  for  a  straight  vortex  line  seg¬ 
ment.  The  wake  model  of  the  inboard  trailing  wake  and  the  shed  wake 
can  therefore  be  constructed  out  of  a  series  of  rectangular,  planar 
vortex  sheets  just  as  the  wake  model  for  the  concentrated  tip  vortex 
is  made  »p  of  straight  vortex  line  segments  (see  Reference  1). 
Breaking  the  wake  up  into  planar  rectangles  adds  artificial  edges 
and  their  associated  singularity.  Since  the  exact  choice  of  the 
boundaries  of  the  rectangles  is  arbitrary,  they  can  be  chosen  to 
avoid  the  singularity  at  the  edges . 

Under  contract  N00019-68-C-0150,  criteria  were  developed  for 
avoiding  the  edge  singularity.  In  addition,  the  subroutine  used  to 
compute  the  Induced  velocity  generated  by  a  planar,  rectangular 
vortex  sheet  was  Improved  so  that  it  now  takes  2.5  to  3  times  as 
long  as  the  corresponding  computation  for  a  straight  vortex  line 
segment  compared  to  the  5  to  6  times  as  long  required  by  the 
original  subroutine.  A  preliminary  trial  of  the  vortex  sheet  wake 
model  was  made  using  the  new  wake  geometry  program.  The  results  were 
promising  but  the  computation  required  almost  3  times  as  long  as  a 
similar  computation  using  the  vortex  line  model.  Further  develop¬ 
ment  is  needed  to  reduce  the  computation  required  and  to  improve 
convergence.  In  addition,  a  better  experimental  data  base  is  needed 
to  properly  evaluate  the  results .  Work  on  these  is  proceeding  under 
contract  N00019-69-C-Q219  and  will  be  reported  in  Reference  2  along 
with  a  detailed  description  of  the  computation  of  induced  velocities 
and  the  criteria  for  u raiding  the  edge  singularity . 
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